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Abstract 

Motivated by the large mixing angle solutions for the atmospheric and solar 
neutrino anomalies, flavor changing neutral current processes and lepton fla- 
vor violating processes as well as the muon anomalous magnetic moment are 
analyzed in the framework of SU(5) SUSY GUT with right-handed neutrino. 
In order to explain realistic mass relations for quarks and leptons, we take 
into account effects of higher dimensional operators above the GUT scale. It 
is shown that the supersymmetric (SUSY) contributions to the CP violation 
parameter in — mixing, ex, the — > 67 branching ratio, and the muon 
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anomalous magnetic moment become large in a wide range of parameter space. 
We also investigate correlations among these quantities. Within the current 
experimental bound of B{n 67), large SUSY contributions are possible 
either in the muon anomalous magnetic moment or in ex- In the former case, 
the favorable value of the recent muon anomalous magnetic moment measure- 
ment at the BNL E821 experiment can be accommodated. In the latter case, 
the allowed region of the Kobayashi-Maskawa phase can be different from the 
prediction within the Standard Model (SM) and therefore the measurements 
of the CP asymmetry B ^ J/ipKs mode and Auib,^ could discriminate 
this case from the SM. We also show that the r ^ /i 7 branching ratio can be 
close to the current experimental upper bound and the mixing induced CP 
asymmetry of the radiative B decay can be enhanced in the case where the 
neutrino parameters correspond to the Mikheyev-Smirnov-Wolfenstein small 
mixing angle solution. 



I. INTRODUCTION 

In order to explore physics beyond the Standard Model (SM), indirect searches play an 
important role complimentary to direct searches of new particles at high energy frontiers. 
The indirect searches include flavor changing neutral current (FCNC) processes, lepton 
flavor violation (LFV) and the muon anomalous magnetic moment. In the minimal SM, 
lepton flavor is conserved and FCNC is forbidden at tree level, so that B, K and /i decay 
experiments have supplied severe constraints on models beyond the SM. At the recent BNL 
E821 experiment, it was reported that the muon anomalous magnetic moment had 2.6 a 
deviation from the SM prediction If the deviation is confirmed by improvement in 
both statistics and understanding theoretical uncertainty of the SM prediction, the muon 
anomalous magnetic moment becomes a clear signal of physics beyond the SM. 

Among candidates of the physics beyond the SM, supersymmetry (SUSY) is the most 
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attractive one. Because of the cancellation of quadratic divergence in the renormalization of 
the Higgs field, SUSY models do not have the hierarchy problem of the SM. Furthermore, 
the gauge coupling unification is realized in the SUSY grand unified theories (SUSY GUT) 
based on SU(5) gauge group or its extensions. 

In view of the flavor physics, it is important that the scalar partners of quarks and 
leptons, namely squarks and sleptons, have a new source of flavor mixing. Due to the new 
flavor mixing, LFV and FCNC such as/i— i>e7, 6-^57 and K*^ — K^/B^ — mixing 
could be induced through SUSY loop diagrams. Because these processes receive too large 
contributions for generic flavor mixing in squark and slepton sector, the structure of SUSY 
breaking sector of the Lagrangian is required to have a special form, unless the masses of 
the SUSY particles are beyond multi-TeV region 0. The simplest possibility to avoid this 
problem is that the SUSY breaking mechanism is assumed to be flavor blind. However, 
even in such a case, the squark and slepton mass matrices receive radiative corrections from 
interactions below the scale where the SUSY breaking is originated, and flavor blindness is 
broken [Q,^. In particular, effects of the large top Yukawa coupling constant can not be 
neglected. A number of analyses have been done in the context of minimal supergravity 
(minimal SUGRA) anzats where SUSY breaking parameters are assumed to be flavor blind 
at the Planck scale . It was shown that the flavor mixing is controlled by the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element. As a result the maximal deviation from the 
SM in the CP violating parameter of — mixing, Sk-, and Bet — Bcl/ B^ — Bg mixing is of 
the order of 10 %, while the SUSY contribution to 6 — > s 7 process can be much important 
|0J^. In the GUT scenario, there are additional contributions to FCNC/LFV processes from 
GUT interactions PJ9|-[T2|. As for LFV processes the /i — 67 branching ratio is close to the 
current experimental bound, especially for SO(IO) model |I0| . 

Recent experimental evidences of neutrino oscillation indicate existence of small neutrino 
mass and large flavor mixings in the lepton sector [jl3| . A natural explanation for small neu- 
trino mass is the seesaw mechanism In this mechanism, heavy right-handed neutrinos 



are introduced, and these neutrinos have Majorana mass term and new Yukawa interactions. 



Because the neutrino Yukawa coupling constants can be as large as the top Yukawa coupling 
constant if the right-handed neutrinos are 0(10^^) GeV, radiative corrections from these in- 
teractions contribute to the renormalization of slepton mass matrix above the mass scale of 
the right-handed neutrinos. Within the minimal SUGRA scenario, it was shown that the 
branching ratios of LFV processes becomes large enough to be measured in near-future ex- 
periments [|T5|-[l7[|. Some GUT models which have predictable neutrino mass and mixing are 



already constrained |18[. In the context of SUSY GUT, these new interactions in the lepton 
sector also contribute to the quark sector, because radiative corrections on the squark from 
neutrino interactions can become a new source of quark FCNC processes as well as the LFV 
processes. Recently, these processes are analyzed in the minimal SU(5) SUSY GUT with 



right-handed neutrino and large deviations from the SM are predicted |ll9| , p0| . However, 
in these analyses, simple flavor structure was assumed so that the correct mass relations 
between the down-type quarks and charged leptons in the first and second generations can 
not be realized. 

Very recently, the BNL E821 experiment reported a new result on the muon anomalous 
magnetic moment [|I|]. The measured value of = {g^ — 2)/2 is a^(exp) = 11659202(14)(6) x 
10^1°, which is compared to the SM prediction a^(SM) = 11659159.6(6.7) x 10"^°. It 
was concluded that the theory and experiment had 2.6 cr difference a^(exp) — a^(SM) = 
43(16) X lO^^*'. The deviation can be explained in the context of the SUSY model pl| . 



In contrast to the LFV and FCNC processes which are very sensitive to the origin on the 
flavor mixing at high energy scale, the muon anomalous magnetic moment can provide us 
information on the slepton masses, rather independent of the flavor structure of the slepton 
mass matrices. 

In this paper, we discuss FCNC/LFV processes in the SU(5) supersymmetric grand 
unified theory with right-handed neutrino (SU(5)RN SUSY GUT) taking account of realistic 
mass relations. The seesaw mechanism generates small neutrino masses and large mixing 
angles which incorporate atmospheric and solar neutrino anomalies. In order to reproduce 
realistic mass relations, we introduce a higher dimensional operator including 24 superfield 



which gives contributions to the Yukawa couphng matrices for the down-type quarks and 
the charged leptons in a different manner. Moreover, new degrees of freedom arise in the 
choice of the bases when the MSSM multiplets are embedded in the SU(5) multiplets. We 
show that the main effect of these new mixings is described by two mixing angles which 
parameterize rotations of the bases between the down-type quarks and charged leptons in the 
first and second generations. We perform numerical analysis on FCNC/LFV processes taking 
account of various sources of flavor mixing. We also calculate the muon anomalous magnetic 
moment and investigate the correlations among these quantities. Solving renormalization 
group equations for the Yukawa coupling matrices and the SUSY breaking parameters, the 
flavor mixing in the squark and slepton sectors is evaluated at the electroweak (EW) scale. 
In addition to the muon anomalous magnetic moment, we calculate following FCNC/LFV 
observables: the branching ratios of /j, ^ ej, t —>■ /ij and b —>■ sj, Sk, the mass differences 
in Bd — Bd mixing and Bg — Bg mixing, and the time-dependent CP asymmetries of 5 — > 
J/ '4> Ks and B ^ Mgj where Mg is a CP eigenstate including a strange quark. We find 
that the SUSY contributions to Ek, the — > 67 branching ratio, and the muon anomalous 
magnetic moment become large in a wide range of parameter space. Within the current 
experimental bound on B(/x — > 67), large SUSY contributions are possible either in the 
muon anomalous magnetic moment or in ek- In the former case, the favorable value of 
the recent muon anomalous magnetic moment measurement at the BNL experiment can 
be accommodated. In the latter case, the new contribution Ek modifies the constraint 
for the CKM matrix elements and affects B decay observables because allowed region of 
Ams^/Ams^ and the time-dependent CP asymmetry of the B J/ipKs mode can be 
quite different from that of the SM or MSSM without the new flavor mixing source. We 
also show that B(r fij) and the indirect CP asymmetry of the radiative B decay can 
be large in the case where the neutrino parameters correspond to the small mixing angle 
Mikheyev-Smirnov-Wolfenstein (MSW) solution. We also notice that the branching ratio of 
/X — > 67 can be close to the present experimental upper limit both in the large and small 
mixing angle MSW solutions. 
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The rest of this paper is organized as follows. In Sec. |T|, the SU(5)RN SUSY GUT 
is introduced. The higher dimensional operators are included to incorporate the realistic 
fermion mass relation. Two new mixing angles are defined to parameterize the effect of these 
operators. In Sec. |T|, the minimal SUGRA model is introduced for the SUSY breaking 
sector. Radiative corrections for the SUSY breaking parameters and FCNC/LFV processes 
are discussed qualitatively using approximated formulas. In Sec. |^ the numerical results 
for the muon anomalous magnetic moment and FCNC/LFV processes are presented. Sec. 
is devoted for conclusion and discussions. In Appendices, useful formulas are collected. 

II. SU(5) SUSY GUT WITH RIGHT-HANDED NEUTRINO 

In this section we discuss quark and lepton Yukawa couplings in the SU(5)RN SUSY 
GUT. Before introducing higher dimensional operators, we first discuss the case without 
them. Later, we introduce those operators to accommodate realistic mass relation. Without 
higher dimensional operators, the Yukawa coupling and the Majorana mass term of the 
superpotential for this model are given by 

Wsu(5)RN = ^e,,,rfe(A„),,,(r)'^''(r^r''/f^+ (A,),,(F),(T^r^:ffb 

o 

+{x,),,Tr(F')aH'' + ^{M,),,'irw, (1) 

where T\ and are 10, 5 and 1 representation of SU(5) gauge group, respectively, i, j 
are generation indices and a, b, c, d and e are SU(5) indices, eabcde is the totally antisymmet- 
ric tensor of SU(5) gauge group. H and H are Higgs superfields with 5 and 5 representations. 
In terms of SU(3) xSU(2)ixU(l)y, contains Q*(3,2,|), f7*(3, and T (1,1,1) su- 
perfields. Here the representations for SU(3) and SU(2) groups and the f/(l)y charge are 
indicated in the parentheses. includes ^(3, 1, |) and L^{1, 2, — |), and iV* is a singlet of 
SU(3)xSU(2)LxU(l)y. H consists of Hc{S, 1, 0) and H2{1, 2, i) and H contains Hci^, 1, 0) 
and Hi{l, 2, — i). {Xu)ij, {Xdjij and {Xu)ij are Yukawa coupling matrices and {M^)ij is a Ma- 
jorana mass matrix. In addition to the above formula, we also need a superpotential for 
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Higgs superfields, Wh{H, H,I1) where E% is a 24 representation of SU(5) group. It is as- 
sumed to develop vacuum expectation values as = diag(|, |, |, — |, —1)vg at the GUT 
scale (Mg ^2x 10^^ GeV) and breaks the SU(5) symmetry to SU(3) xSU(2)i,xU(l)y. 

Below the GUT scale, the heavy superfields such as He, He and E are integrated out 
and the superpotential of the MSSM with right-handed neutrino (MSSMRN) is given by 

Wmssmrn = (yu)ijU'Q'H2 + (ydhjD'Q^H^ + {ye)i,^UH^ + ixH^H^ 

+ {y,)ijN'L^H2 + ^{M,)ijWW, (2) 

where Yukawa coupling matrices are related to those of the SU(5)RN as {yu)ij = i^u)ij: 
{yd)ij = ivDij = {^d)ij and {y„)ij = Below the Majorana mass scale (= Mr), the 

singlet fields are also integrated out from the superpotential and a dimension five operator 
is generated as follows: 

AW. = ~{K,)ij{L,H,){LjH,), = (yJ)^,(_L)^'(y^), .. (3) 

After the EW symmetry breaking, this operator induces by the seesaw mechanism the 
following neutrino mass matrix, 

= (Kh{H,y. (4) 
In this model, the naive GUT relation is predicted at the GUT scale, 

iveh = {yd)ji- (5) 

Although this relation gives a reasonable agreement for and m^, it is well known that the 
mass ratio of down-type quarks and charged leptons in the first and the second generations 
can not be explained in this way. One possibility to remedy this defect is to introduce higher 
dimensional operators above the GUT scale because they can give different contributions 
to the Yukawa coupling matrices of down-type quarks and charged leptons after the SU(5) 
symmetry breaking. 

We consider higher dimensional operators including the 24 Higgs superfield up to di- 
mension five terms. Relevant parts of the superpotential is parameterized as follows: 

7 



AW; 



SU(5)RN 



M 



X 



(6) 



where Mx is the cut-off scale which we take as the Planck mass Mp. We also assume that 
the elements of coupling matrices k^, k^, 7?^ and Ky are smaller than 0(1). After SU(5) 
symmetry is broken, they give contributions of the order of ^ = vg/Mx ~ 0.01 to the 
Yukawa coupling constants of the MSSMRN as follows: 



f 1 5 

{yd)ij = (Ad)ij + ^ ^^{l^d)ij - ^{'^d)r 



(7a) 
(7b) 
(7c) 

(7d) 



The naive GUT relation between the lepton and the down-type quark Yukawa coupling 
matrices in Eq. (R) is modified to 



{.ye)ij = {yd)ji + Q^{'^d)ji- 



(8) 



With this small contribution from the higher dimensional operator, realistic mass relations 
between the down-type quarks and charged leptons can be incorporated in the model. In 



the following analysis we take {k,^ 



because they are not 



necessarily required to reproduce the realistic mass relations. 

In the following, we show that new mixing angles are introduced at the GUT scale because 
of ^d- Using SU(5) symmetry, we can rotate the generation indices of superfields in Eq. (^ 
so that the Yukawa coupling constants and the Majorana mass matrix are parameterized as 
follows: 
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(.Hujij 


= {^CKM^uJi Vuki^CKM) j, 


(9a) 


{ydhj 


= iVD)i'ydj, 


(9b) 


iveh 




(9c) 




= yuAVLYj, 


(9d) 






(9e) 



where Vu, Ve, Vd and are unitary matrices and Vckm is the CKM matrix at the GUT 
scale. Uui, Vdi, yei, y^i, and M^^ represent the eigenvalues of the Yukawa coupling matri- 
ces and the Majorana mass matrix. The GUT relation between the two Yukawa coupling 
constants is then given by 

{Vn)/yci,-ye^{viy^ = la^,).,. (10) 

From this formula wc can derive the following approximate relations for the 1-3 and 2-3 (3-1 
and 3-2) elements of the mixing matrices because the Yukawa coupling constants of the first 
and second generations are much smaller than that of the third generation, 

{VDTs^lUOis, (11a) 
yb 

{VeTs ^ (4)^3, (lib) 

{Ve)\ ^ ^UvS^^dVEh, (11c) 

Qyb 

iVof, ^ -|-(VJ4^z.)3., (lid) 

for i = 1, 2. We can estimate ^/yb and ^/y^ as 

1^1^-1^, + (12) 

yb yr V^rur V2 ^ 

where /9 is a vacuum angle of two Higgs vacuum expectation values (tan/? = {H2)/{H^)). 

Assuming the condition {Kdjij ^ 0(1), we can conclude that the magnitude of these elements 

are constrained to be smaller than (tan/3)"^ because lower tan/3 region is excluded from 

Higgs boson search. On the other hand, 1-2 (2-1) element is not constrained from such a 

consideration. Motivated by this observation we assume the following form: 



{Ve 



e*(-^«+'^) 

e''^^ cos 9e -e-'/^^ sin 9e 
e'f^^sineE e-'^^cosOE 



(13a) 



V 



e*(-^^+'^) 



(13b) 



The antisymmetric part of the Yukawa matrix for the up-type quarks is also written by the 
coefficients of the dimension five operator as follows: 



j ■ 5 -j. 



(14) 



Because we set (k. 



for simplicity, (Vfy)* • becomes e^'^'^^S^^ in our analysis. 



The neutrino Yukawa coupling matrix and the Majorana mass matrix are constrained 
from the oscillation solutions of the atmospheric and solar neutrino anomalies. In the basis 
where the charged lepton mass matrix is diagonal, the neutrino mass matrix is written as 
follows: 



{fnu)ij — ( Kins ) i'^^^i^fcl Kins 



(15) 

where Vmns is the Maki-Nakagawa-Sakata (MNS) matrix At the Majorana mass scale 
Eqs. (S) and (W) are solved as follows: 



m),'y.,iVE)', = ^v^(0J)>'v/^(4j,s) 



(16) 



where O,^ is a complex orthogonal matrix which can not be determined from the low energy 
experiments. Although we neglect a running effect of the neutrino mass matrix between the 
low energy scale and the GUT scale in Eqs. (^) and (|16D, later we fully take account of this 
effect in the numerical calculation in Sec. ITVl. 
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III. MINIMAL SUGRA AND THE MUON ANOMALOUS MAGNETIC 
MOMENT AND THE FCNC/LFV PROCESSES 



In Subsection |III A| we first discuss the flavor mixing of squark and slepton mass matrices 
induced by tlie radiative correction due to tlie Yukawa coupling constants. In order to explain 
qualitative features we show the one-loop logarithmic terms for SUSY breaking parameters. 
In the numerical calculation in Sec. however, we use the full renormalization group 
equation (RGE) and solve them numerically. In Subsection [III B| , we give a brief description 
on the SUSY contribution to the muon anomalous magnetic moment and various FCNC 
and LFV processes. 



A. Minimal SUGRA and radiative corrections to the SUSY breaking parameters 

The soft SUSY breaking terms of the MSSM are given by 

Aoft = -{mlr.QW - {m'uVU-U, - {ml)/mD, 

+ {y,),,E'*UH, + fiBH,H2 + H.c.} 

+^MiBB + ^M2WW + ^MsGG, (17) 

where Q*, W*, D^*, U and E"^* are scalar components of Q*, , -D*, U and respectively. 
We use the same symbols as superflelds for scalar components of the Higgs supermultiplets. 
B, W and G are U(l), SU(2) and SU(3) gauginos, respectively. At the GUT scale, these 
soft SUSY breaking parameters are determined by the following SUSY breaking terms of 
the SU(5)RN SUSY GUT, 
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4 

+^M,g',G,, (18) 

where T', F and A?" are scalar components of T', and 7V* and represents SU(5) 
gaugino. We assume the minimal supergravity scenario for the origin of SUSY breaking 
and set the following boundary conditions for the SUSY breaking parameters at the Planck 
scale, 

{mlYj = {m^y, = {ml,yj = r^S], (19a) 
(X)ij = moAo{X)ij, (A = A„, A^, A^), (19b) 
{K,)ij ^ mo{Ao + AAo){K,)ij, {n ^ , na.'Kd, Hu) , (19c) 
M5 = Mo. (19d) 

If we ignore radiative corrections from the gauge and Yukawa coupling constants and assume 
A^o — 0, the soft SUSY breaking terms are given by 

{mlTj = {mlYj = (ml),^ = {ml)/ = (ml)/ = rr^Si, (20a) 
{y)ij = moAo{y)ij, {y = ya, ye). (20b) 

Then LFV processes are forbidden and SUSY contributions to FCNC processes are sup- 
pressed. We consider AAq ^ case later. 

Radiative corrections between the Planck scale and the EW scale modify the above 
structure of the soft SUSY breaking terms. In particular, the corrections from the Yukawa 
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coupling constants associated with the colored Higgs supermultiplets and right-handed neu- 
trino supermultiplets are important because they have different flavor structure from the 
Yukawa coupling constants of the MSSM. 

Let us estimate these corrections using approximate formulas only considering logarith- 
mic terms to see qualitative features of FCNC/LFV processes in the model. The Yukawa 
couplings including colored Higgs supermultiplets are parameterized as follows: 

Wc = -{ycRhHcU'E' - \{ycL\,HcQ'Q^ 

-{VcR^.TlcD'U' - ivcLhTlcLV' + {ycN)^,HcWW . (21) 

It is convenient to work in the basis where the down-type quark and charged lepton mass 
matrices are diagonal, 

{.yu)ij = yuiiVcKuYj, (22a) 
iVdh = Vdi^i, (22b) 
iVeh = Ve^S], (22c) 
{y.h = y^^iVLyy (22d) 

In this basis, Yukawa coupling matrices in Eq. (^1]) are given by 



iycR)ij = yuiiVcKuVEYj, (23a) 

(ycLh = liiV^KM^'e'^^^^yu, + e^^-^yuAVcKMY,}, (23b) 

{ycRh = yd.e'''"^5}, (23c) 

(ycLh = ye^{v^yp (23d) 

{ycNh=y.AVLVDyj. (23e) 

The radiative corrections to squark and slepton mass matrices from these Yukawa coupling 
constants are approximated as follows: 



Ami ~ ~2 [yiyu + ^yl^ycL + yhd + y^cLyci) (3 + l^on^^o^G 

-2 (yiy^ + y\yd) (3 + \Ao\'')mltw, (24a) 
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Ami ~ -2 {^yuVi + VcrvIr + '^y^Ryhn) (3 + \Ao\'')mltG 

-Ay^yli3+\Ao\^)mltw, (24b) 
Am| ^ -2 (2yeyl + 32/^^1/^^) (3 + \Ao\^)mltG 



-Ay,yl{3+\Ao\^)mltw, (24c) 
Am^ ^ -2 (2yd?/^ + 2ycRy\j^ + ycNVcN) (3 + lAoHm^tG 

-42/,yi(3 + |Aonm2%, (24d) 

Ami ~ -2 {ylve + ^yhiV^L + yiy-) (3 + l^n^o^G 

-2ylye(3 + \Ao\^)mltw - 2yly,{3 + \Ao\^)mltR, (24e) 

where we only take account of logarithmic terms so that tc = ln(|^), tR = -p^p- ln(|^) 
and tw = ln( ^^^^ ). Msusy is a characteristic mass scale of the SUSY particles and 
identified to the EW scale. The off-diagonal elements of the above formulas are sources of 
the LFV and FCNC processes. Keeping only possible large Yukawa coupling constants, the 
off-diagonal elements of the mass matrices are approximated as follows: 

iml)\ ^ -2(Vc^km)%Z/*'(V^ckm)V3 + MmHtc + %) 



:ViKM)\yt5',+6\yt\VcKMr, 
+ (^CKM)'32/t'(^CKM)',} (3 + \Ao\')mltG, (2 ^ j), (25a) 

(ml)/ ^ -6(I^JK:km). V(V^cKM^^)3'(3 + \Ao\')mltG, i^ ^ j), (25b) 



m 



2 ~ _oa/TT/T^ k 
D)i 



-2{V^Vl),'y,,\VlV*r,)A^ + \M')<tG, [i^j], (25c) 



mi)^ ^ -2iyl)\y,,\VLf,{3 + \A,\')mi{tG + t.j), (z ^ j)- (25d) 



m 



q)*j corresponds to the flavor mixing due to the large top Yukawa coupling constant which 
already exists within the MSSM based on the minimal SUGRA. (m|;)/ receives radiative 
correction from the up-type Yukawa coupling constant between the Planck scale and the 
GUT scale. This is a well-known mechanism to induce LFV processes in the SUSY GUT 
Q. We notice that the following important features. 

• There are flavor mixings in the right-handed down-type squark and the left-handed 
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slepton sectors. These mixings are absent in the minimal SUGRA model without 
right-handed neutrino supermultiplet. 

• Because Vl can be related to the MNS matrix, large mixing is possible. 

• The main effect of the higher dimensional operator under assumption of Eq. (^) is 
only rotating the basis of light fermions between dji and sr, cr and /i/j. For each of 
these mixings the rotation is described by one parameter 6'^) in or 6'^; in Vg;. 

In the above discussion, we only considered the radiative correction to squarks and 
slepton mass matrices, however, the trilinear scalar coupling constants, yd and ye also 
receive corrections as follows: 



{yu)ij ~ moAuyuiiVcKuY j 

uiqAq (Amg,)>y„fc(ycKM)^j + yu,{ycKM)\{^rnlf^ 
3 + |Ao|2 ml 

{yd)ij ~ moAaydiS] 

moAo {AmDi^yaj + ydii^mlY ■ 



where Au, Ad, Ae and A A are given by 



(26a) 



3 + |Ao|2 ml 
-^moAA {yd,6] - {VS),'yeM)',} , (26b) 

{ye)ij ^ moAeyeiS] 

moAo {AmDi^yej + ye^AmlYj 



—moAA {ye,6i - {V^),'ydM)'j} , (26c) 



Au^Ao- —gl—tc - —gl—tw + Xylite + tw) + Y.y-l(tG + tn), (27a) 

TTig iO m^Q -^-^ 

168 oMo 88 

Ad^Ao- —gl^tc - —gl—tw. (27b) 
5 mo 5 mo 

168 oMo 48 ^Mo 
Ae^A^- —gl^tc - —gl—tw, (27c) 
5 mo 5 mo 

AA ^ -2Qgl^^tG. (27d) 
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The second terms in the right-hand side of Eqs. (^) are induced by the radiative corrections 
from the Yukawa interactions. The third terms in Eqs. ( |26 b| ) and ( p6c| ) come from the higher 
dimensional operators. They break proportionahty between the trihnear scalar coupling 
matrix and the corresponding Yukawa coupling matrix and generate flavor mixings in the 
left-right mixing mass matrices of squarks and sleptons after the EW symmetry breaking. 
If AAq 7^ in Eq. (|19cj ), we have an extra contribution, AAq to Eq. (p7d| ) and there are 
corrections of order mQAoAAo^tc to Eqs. (^) and corrections of order moAAoC,tx {X = 
G, R, W) to Eqs. (pGf) . Notice that even if we assume AAq = at the Planck scale A A is 
induced by the gauge interaction as shown in Eq. (|27d| ) because the renormalization of 
and Kd are different due to the wave-function renormalization of S. We therefore expect 
analysis with AAq = gives us a qualitative feature for general cases. 



B. The muon anomalous magnetic moment, FCNC and LFV processes 

Let us discuss the muon anomalous magnetic moment, FCNC and LFV processes in the 
model according to the approximations of the previous section. In the following we make 
a simplification in the neutrino sector to estimate possible deviations from the SM in the 
FCNC/LFV processes. We assume the Majorana masses of the right-handed neutrinos are 
universal at the Majorana mass scale Mr as {M,y)ij = dijMn. We also neglect any CP 
violating phase in the model except for the Kobayashi-Maskawa phase and assume that Vn, 
Vl, Ou and Vum are real matrices and aD,E = Pd,e = 1d,e = 5 = in Eq. (|T3|) and (pui = 
in Eq. If we include these phases, new contributions to the electron and neutron electric 

dipole moments (EDMs) are induced so that we have to take into account constraints to 
SUSY parameters from the upper bound of the EDMs. In the numerical calculation we 
evaluate these EDMs and check that these constraints are satisfied in the case that new 
CP phases are set to vanish. With the above simplification, the mixing matrix Vl and the 
neutrino Yukawa coupling can be related to the low energy observables according to Eq. (|16|) , 

VL = Vi^s, yu, = ^MRm,J{H2). (28) 
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We parameterize the MNS matrix assuming maximal mixing for the atmospheric neutrino 
oscillation as follows: 



MNS 



COS 6^ 



sin 

cos esni 







(29) 



V2 V2 V2 

sin gsun COS gsun 1 

V ^2 V2 V2 

where 6'sun is the mixing angle for the solar neutrino oscillation. We assume the 1-3 element 
of the MNS matrix is zero in our analysis because it is known to be small from the result 
of the CHOOZ experiment We will comment on the nonzero case later. We assume 
the hierarchical pattern of neutrino mass, namely mj,i < mi,2 <^ mi,^. From the following 
relation, 



2 A 2 I 2 



m 



(30) 



where Arn^^^ and Aml^^ are the mass differences of solar and atmospheric neutrino oscil- 
lation, m^^ and m^2 ^i-re determined once we fix m^i. Then using Eq. (pH]) we can calculate 
yui for a fixed value of Mr. 



1. The muon anomalous magnetic moment 
We consider the SUSY contribution to the muon anomalous magnetic moment 



2T|J2112^ . The muon anomalous magnetic moment is defined by the following effective 



Lagrangian: 



where e is the positron charge, is the muon mass , F^p is the electromagnetic field 
tensor and = 2[7a,7/3]/2. The SUSY contribution to (= a^^^"^) is obtained from the 
flavor diagonal parts of photon-penguin diagrams including smuon-neutralino and sneutrino- 
chargino. a^^^^ depends on the slepton mass and the neutralino/chargino mass and mixing, 
but it is rather insensitive to the flavor mixing of the slepton sector. Therefore we expect 
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almost the same as a result in the minimal SUGRA model. In the MSSM based 
on the minimal SUGRA, it is known that the main contribution comes from the sneutrino- 
chargino diagram which contains a component proportional to /itan/3 [^. Then a^'"^^'^ 
preferred by the recent results from BNL E821 experiment is achieved in the large tan/? 
region of the parameter space. In this region the sign of a^^^"^ is correlated to the branching 
ratio of 6 — > s 7 through the sign of the Higgsino mass parameter /i so that positive 
when 6 — > s 7 is suppressed and negative when 6 — > s 7 is enhanced. 



We consider LFV decays of charged leptons. Radiative decays of charged leptons oc- 
cur through photon-penguin diagrams including sleptons, neutralinos and charginos. The 
effective Lagrangian for these processes is described as follows: 

^LFv ^ {m,,A'^(l^,a'^''lLj)F^, + me,A'i{&''lnj)F^,} + H.c. {t > j), (32) 

where Gp is the Fermi constant and i,j denote generation indices. corresponds to the 
amplitude for If I'jlR and /C'l for I'l The branching ratios are calculated from 

these amplitudes as B(/+ ^ /+7) = ?,M'k'^{\A'1\'^ + \A'l\^). 

For /i — 6 7, it is known that if both left-handed and right-handed sectors have flavor 
mixing, there are the diagrams which have an enhancement factor nij. as shown in Fig. |l| 
10| , P^JT7| . In our model we find that diagrams corresponding to Figs. ^ and ^ also give 



large contributions. The flavor mixing in the left-right mixing term in Fig. |^ is induced by 
renormalization between the Planck and GUT scale as shown in Eq. (|26c|) . Approximate 
formulas for A^jf and A^jf from these contributions are given by 



Ar ~ -^sin20,un(l/i.2^ - VuiVt 



Vckm)\ sin Oe + (Vckm)^2 cos Be} 
— cos9Esin9D a^, (33a) 
~ -{yu3^ - yu2')yt {{Vckm)\cos9e - {VcKM)\sm9E} —4 + sm9ECOs9D <, (33b) 
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where we explicitly show ^sun, Uui, Ut, Oe and 6d dependence, a^, a'^ and a" are functions of 
the slepton masses and the chargino and neutralino masses and mixings. These contributions 
correspond to Fig. |1|, Fig. ^ and Fig. ^ respectively. The explicit forms of the functions are 
given in Appendix |^. Because (Vckm)^2 ^ (K:km)^i, the mixing angle 9e can enhance the 
amplitude Aj^ compared to the case 9e = 0- 

The ratio of the magnitudes of A^^ and A"^ can be measured by the P-odd asymmetry 
of /i ^ 67 process, ^67) ||26|. With the help of initial muon polarization, we define 
A{fj, — i> 67) as follows: 

"^^^^J^f = ^B(/.+ e+7) {1 + ^(/^+ - e+7)i^cos^} , (34a) 
I 42112 _ I 42112 

Mf^^^e^l)= \Al , Al ^ (34b) 



where P is the polarization of initial and 6 is the angle between the polarization and the 
momentum of the decay positron. For r — > yU7, a similar P-odd asymmetry can be measured 
in the e^e^ r+r^ process using spin correlation of the r pair 



3. h^s-i 

The A5 = 1 FCNC effective Lagrangian for the radiative B decay is written as follows: 

£^^=1 = --^ {C',(Ma^''hL)F,, + C^(JZa^'%R)F,,} + H.c. (35) 

In the SM case, the process occurs through photon penguin diagrams which exchange a W 
boson as Fig. ^ and C7 is suppressed by a factor of rris/mb compared to C7. In the minimal 
SUGRA model without the right-handed neutrino supermultiplets, the flavor mixing in the 
squark mass matrices only appears in that of left-handed squark and the same argument can 
be applied. In the present model, however, a gluino exchanging diagram (Fig. |^) can give 
a large contribution to Cj because of the new flavor mixing in the right-handed down-type 
squarks {m'jj)2^- 
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If C7 has a similar magnitude as C7, the time-dependent CP asymmetry oi B ^ Ms 7 
may be observed where Mg is a CP eigenstate which includes a strange quark such as Ki 
(— > Ksp^) or K* {-^ KsTi^) |T9|,|28[|. The asymmetry is defined as follows: 



where T{t) {T{t)) is the decay width of B^{t) ^ M^-f (5°(t) M37). is +1 if M, is a 
CP even state and —1 if Mg is a CP odd state. 6b is the phase of B^ — B^ mixing amplitude 
Mi2{Ba) which defined below in Eq. (^Ol). In the SM case, this asymmetry is only a few 
percent, however, it may be considerably enhanced by the new SUSY contribution to Cj. 



4. £k 

— mixing is described by the AS* = 2 FCNC effective Lagrangian. The general 
form is given by 

1 1 

+ 9RL(dlsRa)(dRSLp) + 9liL(dlsRp)(d1^SLS^ + H.C., (37) 

where a and (3 denote color indices. The explicit forms of effective coupling constants in the 

above formula are given in Appendix 0. The CP violation parameter in — mixing, 

Ek is calculated from the above formula as follows: 

ef- Im{M,,(ir)} {K'\C^^='\K') 
= 71 Am^ ' ^''^^^ = 2^ • ^^^^ 

In the case of the SM, the process occurs through a box diagram in which two W bosons are 

exchanged between the down-type quarks (Fig. |^) so that it is dominated by the gfj^-term. 

However, we have a new flavor mining, {mj^)^^ in the right-handed down- type squark sector. 

We can draw gluino exchanging diagrams as Fig. which include the CP violating phase 
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of CKM matrix in (?71q)^2 of ^^e squark lines and large flavor mixing in {fnj^)^^ 

on the other. These diagrams contribute to the coupling constants g^j^ and gfn^' ■ In the 
actual numerical calculation we first derive the effective Lagrangian at the energy scale 
MgusY- According to the reference [^], we include QCD corrections and derive the effective 
Lagrangian at the hadronic scale. The matrix elements for the dominant operators are 
parameterized as follows: 

{K'\(dh^SLa)(dll,SLp)\K') = (39a) 
{K\T^s:^^){dlsn,)\K') = 1 (^^)' </|(5^)|^, (39b) 

{K\dls,,)(disna)m = I (^^y myUBKruL, (39c) 

\ms + TTld/ 

where B^, (-Bx)'rl and (-Bx)'rl' bag parameters calculated by the lattice QCD method. 
Because there is a large enhancement factor of order (m/^/ms)^ in the matrix elements and 
large mixngs originate from the MNS matrix, SUSY contribution to Ek is expected to be 
large in this model. 

5. Bd - Bd/Bs - Bs mixing 

The Ai? = 2 effective Lagrangians for Ed — Bd and Bs — Bg mixings are parameterized 
in the same manner as — mixing. For Bd — Bd mixing, it is obtained by replacing the 
strange quark with the bottom quark in Eq. (pTj). For Bs — Bs mixing, we further replace 
the down quark by the strange quark. The mass difference of Bd — Bd mixing, Aitlbj^ is 
calculated from the effective Lagrangian as follows: 

ArriB, = 2\MuiBd)\, M^iBd) = -^-^ (40) 

2ms^ 

In the SM case, this process occurs through the W boson exchanging diagram and the g^- 
term gives a dominant contribution. In the present case, there are new diagrams as shown in 
Fig. ^ which contain the new flavor mixing in (m|))i^ on one of the down- type squark lines 
or on the both of them. The former contributes to gf^j^ and gf^j^' and the latter to g^. Unlike 
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— mixing, the scalar-scalar matrix elements do not have an enhancement factor for 
Bd — Bfi mixing case because mB^/mi, is 0(1). Similar argument holds for — Bg mixing. 
In the numerical calculation we use next leading order QCD correction for g\ and and 
leading order QCD formulas for other contributions ||2^ and bag parameters are calculated 
by the lattice QCD. Numerical values are shown later. 

IV. RESULTS OF NUMERICAL CALCULATIONS 

In this section we present our numerical results on the FCNC/LFV processes in the 
SU(5)RN SUSY GUT. 

In the present analysis we assume that the SUSY breaking terms have the minimal 
SUGRA type boundary condition at the Planck scale and that the Kahler potential is flat. 
Adopting the simplifications discussed in the previous section, we have the following input 
parameters. 

• Parameters at the Planck scale: the universal scalar mass mo, the universal gaugino 
mass Mo, and the universal coefficient for the scalar couplings Aq. 

• Parameters at the GUT scale: mixing angles 6d and 6e- 

• Parameter at the right-handed neutrino mass scale: Majorana mass of the right-handed 
neutrino Mr, which is also used as the matching scale. 

• Parameters at the EW scale: quark, lepton and neutrino masses, mixing matrices 
VcKM and Vmns, tan/5 and the sign of the Higgsino mass parameter in Eq. (g). 

Throughout the following calculation, we fix some of the parameters as shown in Table |. 
We consider two cases for the neutrino parameters, corresponding to the large mixing angle 
(LMA) and the small mixing angle (SMA) MSW solutions of the solar neutrino anomaly. 
The parameters we used in the neutrino sector for each case are given in Table For Mr 
and tan/?, we take several cases to see the dependences (see Table fT|). SUSY breaking 
parameters mo, Mq and Aq are varied and AAo is fixed to zero. 

22 



With these parameters, we solve the RGEs of the mass parameters and the couphng 
constants between the Planck and the EW scale taking all the flavor mixings into account. 
Detail of our method is explained in Appendix The magnitude of fi is determined by the 
radiative EW symmetry breaking condition, in which the minimum of the one-loop effective 
potential for the Higgs fields is evaluated. Then we obtain all the masses and mixings of the 
SUSY particles at the EW scale and calculate the FCNC/LFV observables as functions of 
above parameters. We calculate the following quantities: 

• The SUSY contribution to the muon anomalous magnetic moment a^^^^; 

• Branching ratios of6-^s7, >e7 and r /i 7; 

• P-odd asymmetry of /i — >• 67; 

• — mass splittings ^fnB^ and Aitlb/, 

• CP violation parameter ek', 

• Time-dependent CP asymmetries of B ^ Ms'y and B J/ipKs- 

In order to find the allowed region in the parameter space, we impose the constraints from 
the experimental results of the direct searches of SUSY particles ||3^ and Higgs bosons ||31 



and the measurements of B(6 — > S7) |3^. Also it turns out that, in some parameter region, 
the branching ratio of /x — e 7 exceeds the present upper limit and hence this process already 
gives an important constraint on the parameter space. We discuss the constraints from the 
measured values of Ek and Am^^ and the lower bound of Am^^ later, since it depends on 



the CKM parameters, namely \Vub\ and ^13 p3[] 



A. 9e and Od dependence of /i ^ 67 and 

Let us first discuss the 9e and 9d dependence of the ^ 67 decay and Ek- 
As given in Eq. (P5D, the decay amplitudes Aj^ and depend on 9e and On differently, 
so that both of the branching ratio and P-odd asymmetry are affected. In Fig. ^ we show 
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B(|U — > 67) and A{fi — 67) as functions of 9e and 9d- The shaded regions are excluded 



by the upper bound of B(yU — > 67) [Q. Here we take the LMA case for the neutrino 
parameters, = 4 x 10^^ GeV, tan/3 = 20, /i > 0. SUSY breaking parameters are also 
fixed as Mq = 300 GeV, Aq = and mo = 0, 300, 600, 900 GeV. For the fixed 60 = case 
((a) and (b)), we can see that the amplitude A^^ is enhanced for a nonvanishing 6e and 
relatively small mo- In the parameter region 6e ~ 90°, B(yU -^67) becomes larger than that 
for 9e ^ and A(fi -^67) approaches to +1, reflecting that Aj^ is enhanced and dominates 
over A"^. For 6e = case ((c) and (d)), A"^ dominates in the most of the range of 6'^ and 
hence A{fi 67) is close to —1. In some special case, 9d = —30° and uiq = 300 GeV 
for example, a cancellation among contributions to A"^ occurs and the branching ratio is 
suppressed. In such the P-odd asymmetry approaches to +1. 

Fig. |1^ shows the 9d dependence of Ek for the same parameter set as Fig. |](c) and (d). 
This dependence comes from (7^^ and gf^^' in Eq. ( p7| ) since 6e> directly affects the mixing 
between the right-handed down-type squarks of the first and the second generations. We 
have checked that 9e dependence is negligible for Ek- 

Hereafter we fix 6'^; as = and in most cases we also fix O^) = 0. 

B. a?,^^^, FCNC and LFV observables for different sets of the neutrino and the 

SUSY parameters 



In Fig. |TT| we show contour plots of the SUSY contribution to the muon anomalous 
magnetic moment a^^^"^, branching ratios of /i ^ 67 and r ^ /i7, and the deviations 
from the SM values of Ek, ^f^B^ and ArriB^- The input parameters used in each figure 
are given in Table fTT[ We also fix the sign of /i as /i > in these figures, /i < region 
is disfavored because the SUSY contributions to 6 — > 57 decay amplitude interferes with 
the SM contribution constructively, so that the branching ratio becomes too large in a large 
portion of the parameter space. Note that the constraints from Ek, ^^Ba Am^^ are 
not imposed in Fig. [TT| since these constraints depend on 613 and \Vub\- We show our result 
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for each of these observables by taking a ratio to the corresponding SM value, expecting that 
the most of the dependences on the CKM parameters cancel. In fact we have checked that 
the plots do not change when a different value of 613 is used. The dependences of a^^^^, 
B(yU —>■ e'-f) and B(r ^ fi'-f) on the CKM parameters are also small. 

In Fig. |TT|(a) we take the LMA case for the neutrino masses and mixing, = 4 x 
10^^ GeV, Oe = Od = 0, tan/5 = 20 and = as a reference point. Shaded regions 
are experimentally excluded region. The constraints mainly come from the LEP II Higgs 
boson search and the upper bound on B(/i — > 67). We see that there is a parameter region 
with 20 X 10^1° < a^usY < 60 x 10~^°, which is favored by the E821 result and in that 
region B(/i — 67) becomes larger than 10"^^. In the allowed parameter region within the 
plotted range Mq < 1 TeV and mo < 4 TeV, B(/i ^67) varies 0(10"^^) to O{10-^^) and 
B(r ^ fi'y) varies 0(10^^^) to 0(10^^). Both branching ratios depend similarly on Mq and 
mo. Also we see that the deviation of Ek from the SM value is about ten percent at most 
and the deviations of ^itlb^ and Am^^ are small. 

In Fig. 0(b) plots for Aq = 2 are given. In this case the stop mass squared becomes 
negative in some parameter region, which is shown in the figure. The excluded region by 
the B(/i — i> 67) constraint is enlarged, due to the enhancement of the branching ratio by 
the change in the left-right mixing in the slepton mass matrices. The deviation of Ek is also 
enhanced and there is an allowed parameter region where Ek is enhanced by more than 25 
percent. It is noticeable that the allowed parameter region with a large enhancement of Ek 
is different from the E821-favored region. The region which corresponds to both a favorable 
a^^^^ and a large enhancement of Ek is excluded by other constraint, such as B(/i — 67). 

Comparing Fig. |ll](a) and (c), we can see the dependence on 9d- Since 9^ affects 
the mixing between the first and the second generations, the difference appears mainly for 
B(/i — > 67) and Ek- For a nonvanishing 6^ = 45° (Fig. 0(c)), B(yU —>■ e'-f) is enhanced for 
mo ^ 700 GeV and the excluded region is enlarged. Also the SUSY contribution to Ek can 
be larger than the SM contribution in a part of the allowed parameter region as shown in 
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Fig. |Ty. a^^^ does not depend on so that the parameter region with a^^^ ^ 10 x 10 



10 



is excluded in the case (c). The behavior of B(r ^ fi'j) and ArriB^ are unchanged also. 
Although ArriB^ depends on 6^, the deviation is quite small in either case. 

The plots for tan/5 = 5 are given in Fig. |ll|(d). Since both B(/i — > 67) and B(r — > yU7) 
are proportional to tan^/5, possible values are suppressed as B{fi ^ 67) < 10^^^ and 
B(r — /U7) ^ 10^^°. a^^^"^ is proportional to tan/5 and is also suppressed. The excluded 
region is larger than the tan /5 = 20 case because the constraint from the Higgs mass bound 
is stronger for a smaller tan/3. The plots for ek, ^f^B^ and Am^^ are the same as those in 
the case (a) except that the excluded region is enlarged. 

Fig. |Tl|(e) shows the case with a larger Mr = 4 x 10^^ GeV. Since the magnitude of the 
neutrino Yukawa coupling constants are proportional to y/Mji, the flavor mixings in m|, 
and m| are enhanced for a larger M^. As a result we see that B(yU — 67), B(r fi'y) 
and Ek are significantly enhanced in this case, compared to the = 4 x 10^^ GeV case 
(a). Although the excluded region due to the constraint from B(/i — 67) is enlarged, there 
is still an allowed parameter region where Ek is enhanced more than fifty percent of the 
SM value. Also B(r ^ ^^7) can be close to O(10~^). In the allowed parameter region, the 
deviations of Am^^ and Am^^ are small, unaffected by the change of Mr and the 

E821-favored region is excluded by the B(/i —^67) constraint. 

Fig. |Tl|(f) shows the SMA case. Other parameters are taken to be the same as those 
in the case (a). In this case, the mixing between the first and the second generations is 
suppressed compared to the LMA case. Consequently B(/i ^ 67) is at most 0(10^^^) in 
the allowed region and the deviation of Ek is smaller, a^^^'*', B(r Ai7), Am^^ and Aitlb^ 
look the same as those in the case (a). 

In all the above cases the deviation of the — mixing from the SM value is small. 
Let us now show an example with a large enhancement of the B^ — Bg mixing in Fig. |ll](g). 
We see that Aitlb^ differs from the SM value by more than 50 percent in a parameter region 
Trio ~ 700 GeV and Mq ^ 200 GeV. In the same region ek is also enhanced by a similar 
amount. Note that this enhancement comes from the mixing in the right-handed down- type 
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squarks induced by the neutrino Yukawa coupling. In such a case the time-dependent CP 



asymmetry oi B ^ Mg'y decay is also enhanced. Fig. |T2| shows Acp{B M^^y) with 
the same parameter set. We see that this asymmetry can be larger than 25 percent in the 
parameter region where ArriB^ is enhanced. In this case B(r —>■ fi'j) can be close to 10~^. 

In Table |V| we summarize the possible SUSY contributions to the observables given in 
Fig. |ll]. We can see that, except for the case (g), a large deviation from the SM is possible 
only in a^^^ , B(/i — > 67) and ex- 
Let us see the correlation among a^^^^, B(yU ^67) and Ek more closely. Fig. |13| shows 
the correlation between a^^^ and B(/i — > 67), a^^^^ and sk/ {£k)sm, and eK/{£K)sM and 
B(yU —>■ e'-f). Here SUSY breaking parameters mo, Mq and Ao are scanned within the range 
mo, Mq < 3 TeV and —5 < Aq < 5. Other parameters are taken to be the same as those 
in Fig. |ll](a) and (b). In the plot of the correlation between a^^^"^ and B(yU ^ 67), The 
a^^^^ < branch corresponds to /i < and a^^^^ ^ —20 x 10~^° region is excluded by the 
B(6 ^57) constraint. Notice that the parameter region where a^^^^ saturates the E821 
result is different from that with a large £k/{£k)sm- As can be seen in the plot of a^^^ 
and eK/{sK)sM, when ek is enhanced by ~ 50 percent, the magnitude of small. 

C. Allowed region of AtubJ Ams^ and the CP asymmetry of B ^ J/t/j Ks 

Finally, let us discuss the effect of varying the CKM parameters and ^13. Within 

the SM, these parameters are determined by combining the measurements of several observ- 
ables: b uiu semileptonic decays, Ek, Am^^, AnisJ Atjib^ and the time-dependent CP 
asymmetry oi B J/ip Ks decay. However, in the present case we have shown that there 
can be a significant SUSY contributions to these observables, especially for ek- In such a 
case the allowed range of ^13 given by the measured value of Ek is different from the SM 
case, and then this change affects the other observables. As for |\4b/Vc;,|, we have no change 
since the b —>■ uiu decay is dominated by the tree-level SM amplitude. 

We show how this effect will be observed in Fig. IH, where the possible region in the 
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space of Attib J AniBj^ and the time-dependent CP asymmetry oi B ^ J/ip Ks is presented. 
In this figure, we vary iKib/V^fe] and 5i3 within the ranges 0.08 < iKfc/Kbl < 0.1 and 
< 5i3 < 360°. The dotted fines in eacfi plot sfiow tfie SM values of AttlbJ Ams^ and 
Acp{B J/ipKs) for tfie wfiole range of 5i3 and |\4fe/K;fe| = 0.1 (outer line) and 0.08 
(inner line). Tfie sfiaded region is allowed in tfie SM case. We impose tfie constraints from 
tfie measured values of Ek = 2.28 x 10^'^ and AniB^ = 0.482 ps"^ and from tfie lower limit 
of Atub^ > 14.3 ps~^. In tfie calculation of Ek, Aitlb^ and Atub^ we fix tfie bag parameters 
and tfie decay constant of tfie B meson fBa,Bs as given in Table |3^. Wfien we impose 



tfie experimental constraints, we allow ±15 % and ±40 % deviations for ek and Aitlb^, 
respectively, in order to take tfieoretical uncertainties in tfie bag parameters into account. 
Since tfiis uncertainty is expected to be reduced in tfie ratio AitibJ AniB^, we use tfie lower 
limit of tfie ratio AtubJ Aitib^ instead of Aitib^ itself. 

Fig. 0(a) sfiows tfie result witfi tfie same parameter set as Fig. |ll](a) except tfiat Aq is 
scanned witfiin —5 < Aq < ±5 (see Table PTT] ). As sfiown in Fig. |Tl](a) and (b), tfie SUSY 
contributions to B^ — B^ mixing and B^ — Bg mixing are quite small in tfiis case so tfiat tfie 
allowed region lies between tfie dotted lines. Tfie difference of tfie allowed regions from tfie 
SM one comes from tfie fact tfiat tfie SUSY contribution to Ek can be as large as 50 percent 



of tfie SM value, wfiicfi can be seen in Fig. [13 



In Fig. |T|(b) we take Od = 45° as in Fig. |Tl|(c) and Aq is scanned witfiin —5 < Aq < 
±5. In tfiis case tfie enfiancement of Ek is more significant compared to tfie Ob, = case. 
Consequently a region witfi smaller 613 is now allowed and fience a smaller Acp{B — >■ 
J/ipKs) ~ 0.4 is possible, compared to tfie SM value Acp{B — » J/ipKs) ~ 0.7. At 
tfie same time AitlbJ AmB^ can be as large as 60. 

Fig. 0(c) is tfie case corresponding to Fig. |ll](g) witfi —5 < Aq < ±5. In tfiis case 



tfie allowed region can be outside of tfie dotted circles, since a large deviation of Auib^ is 
possible. On tfie otfier fiand, we see tfiat tfie deviation of Acp{B J/ip Ks) from tfie SM 
value is small. 

At present we only fiave lower bound for Auib^ and tfie CP asymmetry of 5 — > J/ip Ks 
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and related modes is not precise enough p5| , |37| . In a few years we expect that the Ams^ will 
be measured at Tevatron and the precision of the CP violating asymmetry will be improved 
to 10 percent level at Belle, BaBar and Tevatron experiments. It is conceivable that the 



deviation shown in Fig. O will be clearly seen in these experiments. 



V. CONCLUSION AND DISCUSSIONS 

In this paper we have studied the FCNC and LEV processes as well as the muon anoma- 
lous magnetic moment in the framework of SU(5) SUSY GUT with right-handed neutrino 
motivated by the large mixing angle solutions for the atmospheric and solar neutrino anoma- 
lies. In order to explain realistic mass relations for quarks and leptons, we have taken into 
account effects of higher dimensional operators above the GUT scale. It has been shown that 
there appear new mixing angles in the right-handed charged leptons and the right-handed 
down-type quarks due to the higher dimensional operators. We have calculated various 
low-energy observables by changing parameters of the model, namely SUSY parameters, 
neutrino parameters (LMA or SMA, and Mfj) and the above new mixing angles. We have 
shown that, within the current experimental bound of B(/i ^67), large SUSY contribu- 
tions are possible either in the muon anomalous magnetic moment or in Ek- The parameter 
regions which have a large correction in one case is different from that in the other case. In 
the former case, the favorable value of the recent result of the BNL E821 experiment can be 
accommodated. In the latter case, the allowed region of the Kobayashi-Maskawa phase can 
be different from the predictions within the SM, and therefore the measurements of the CP 
asymmetry of B ^ J/ipKs mode and Am^^ can discriminate this case from the SM. We 
also show that the t ^ fi'-f branching ratio can be close to the current experimental upper 
bound and the mixing-induced CP asymmetry of the radiative B decay can be enhanced in 
the case where the neutrino parameters correspond to the small mixing angle MSW solution. 

Finally there are several remarks. 

• In this paper we have neglected the constraint from the nucleon decay. If we take the 
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minimal model for the Higgs sector at the GUT scale, it is likely that the nucleon 
decay experiments excludes most of the parameter space even if the squark mass is 



multi-TeV p8[. It is known, however, that there are several ways to suppress the 



nucleon decay without changing the flavor signals discussed here [3S 



For LFV search, /i — > e e e decay and fi — e conversion in a muonic atom are also 
promising experimentally EO] . The rates of these processes have simple relations with 



B(/i ^ 67) if the photonic operator Eq. (|3^ ) gives dominant contribution |[T7|j4T 



B(/i+ ^ e+e+e-) a , mf, 11 



B(/i+^e+7) Svr 
B{fi~N ^ e- N) _ B{A,Z) 
B(/i+ ^ e+7) ~ 428 



2 

0.006, (41a) 



log — I - ^ 

mi 4 



(41b) 



where B{A, Z) represents the rate dependence on the mass number A and the atomic 
number Z of the target nucleus: B{A, Z) ^ 1.1 for ^^Al, B{A, Z) ^ 1.8 for "^^Ti and 
B{A, Z) 1.25 for ^°'^Pb. These relations hold also in our case. 

We also calculated r — 67 branching ratio. In all cases B(r ^ 67) is smaller by two 
or three orders of magnitude than B(r yU7)- 

As shown in Fig. 0(a) and (e), the flavor mixing effect due to the neutrino Yukawa 
coupling is enhanced (suppressed) for a large (small) since the neutrino Yukawa 



coupling constants are proportional to yM^ for given neutrino masses. When we take 
a small value of M/j, such as Mr < 10^° GeV, the contributions of yy in m\ and 
m\) given in Eqs. (|25c|) and (p5d| ) are suppressed and hence the SUSY contribution 



to Ek becomes smaller than ~ 10 percent. Even in this case, however, there are 
contributions to the /i — 67 decay amplitudes independent of the magnitude of y^, 
as shown in Eq. (|33D . The terms proportional to a" dominate the amplitude for 
Od,e = 0(1) and the branching ratio can be as large as the experimental upper bound 
in some parameter region. 
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We have assumed 1-3 element of the MNS matrix to be vanishing. However, present 
experimental upper bound is given as sin^ 29i3 < 0.1 |2^. When a nonvanishing Oi^ is 
introduced, fi —>■ e'-f and Ek are generally enhanced since the loop diagrams including 



the third-generation squarks/sleptons in the internal lines give large contributions |]T7 . 
Consequently, the constraint from the upper bound of B(/i — >■ 67) is significant even 
in the SMA case. In the allowed region, a^^^^ and the SUSY contributions to Am^^ 



and AniBg are smaller than those in the ^13 = case shown in Fig. [Tl|. We see that 
the large deviation of AitlbJ Ams^ outside of the dotted lines given in Fig. [l^(c) 
disappears when we take sin^2^^i3 ^ 0.001, and the corresponding plot looks similar 
to Fig. [T|(b). 

Let us now discuss about the validity of the simplification imposed in the mixing 
matrices Vd and Ve- We have numerically checked that, when we require < 4 

for example, the mixing angles for the second-third and first-third generation mixings 
are restricted to be smaller than ~ 15° for the tan /5 = 20 case. In this case B(yU — > 67) 
varies within the range which is several times larger than those shown in Fig. |^. In 
addition to the mixing angles, CP- violating complex phases of 0(1) can be introduced 
in Vd and Ve- It turns out that the SUSY contributions to ek can be twice as 
large as those given in Fig. 0. These complex phases also contributes to the EDMs 
of the neutron (dn) and the electron {(1^). We calculated EDMs and obtained that 
\dn\ < 10-26 e cm and |4| ^ 10"^^ e cm for mo = 600 GeV, Mq = 300 GeV and Aq = 0. 
Thus the EDMs can be close to the present upper bounds \dn\ < 6.3 x 10"^^ e cm |^ 
and \de\ < 4.0 x lO^^^ g cm |3|. 



As discussed above, if we relax the simple assumptions for the mixing matrices Vd, Ve and 
Vmns, typical patterns of the deviation from the SM can be summarized in the following way. 
(1) B(yU — » 67) can be close to 10"^^ and the deviation in AmBjAmB^-Acp{B J/ipKs) 
plane appears like Fig. |14|(b). In this case, a^^^^ is quite small and the SUSY contribution 
does not saturate the observed discrepancy of a^. (2) compatible with the E821 
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result and B(/i ^67) can be as large as 10~^^. However, no deviation may be seen in ek, 
Bd — and Bg — Bg mixings in this case. From these observations we can conclude that it 
is important to search for new physics effects in the ongoing and near-future experiments, 



namely the BNL muon g — 2 experiment, /i — > 67 and yU — e conversion experiments p^ , ^ 
B physics experiments at 5-factories and Tevatron. Combining results obtained in these 
experiments we may be able to get some insights on interactions at the GUT or right-handed 
neutrino scales. 
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APPENDIX A: RGE AND MATCHING CONDITION AT THE GUT SCALE 

In this appendix we show the detail of our numerical calculation taking account of the 
effects of higher dimensional operators. An outline of the calculation is as follows: 

• We solve the RGEs for the gauge and Yukawa coupling constants between the EW 
scale and the GUT scale. The neutrino Yukawa coupling constants are calculated with 
Eq. (pSD at the Majorana mass scale. 

• At the GUT scale, the coupling constants in the superpotential of the SU(5)RN SUSY 
GUT are determined from the Yukawa coupling constants for quarks and leptons using 
the matching condition explained in Subsection [A2|. Then we solve the RGEs for these 



constants between the GUT scale and the Planck scale. 

At the Planck scale, we set the boundary conditions for the SUSY breaking parameters 
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as Eq. (|T9|) and solve the RGEs for these parameters between the Planck scale to the 
EW scale. 



The RGEs for the MSSM and MSSMRN are given for example in the reference P,p!7[] and 



we show the RGEs for the SU(5)RN SUSY GUT in Subsection In Subsection ^ 



we 



explain the matching condition at the GUT scale taking account of the higher dimensional 
terms in the Kahler potential. 



1. RGE for the SU(5) SUSY GUT with right-handed neutrino 

In this subsection we show one-loop RGEs for the SU(5)RN SUSY GUT. In the derivation 
of RGEs, we only take account of diagrams to which the higher dimensional operators are 
inserted at most one time. In this approximation, the quadratic divergence does not appear 
in the calculation. 

The RGEs for the SU(5) gauge coupling constant and the gaugino mass parameter 
M5 are given by 

(4vr)'M^<75 = hgl (Ala) 
(47r)2M^M5 = 2hglM^, (Alb) 

where M is the renormalization scale. The coefficient of the beta function is given as 65 = —3 
for the minimal field contents. 

The RGEs for the coupling constants in the superpotentials, Eqs. (|I|) and (^) are repre- 
sented as follows: 



(47r 


fM (A„ 


'ij 


= (A«; 




\ + M 






(A2a) 


(47r 




'ij 




kjiQ-p) 


\ + (A.) 




+ (Ad)y0-ff, 


(A2b) 


(47r 




'ij 






)\ + (A.; 






(A2c) 


(47r) 










)\ + {^^t 




)^ + ('^^)..(0H + es), 


(A2d) 



33 



d_ 

\2 : 



d 



where G's are given by 

94 

(eF)^• = 4(AS)^^Al)., + (At)^'^(A.)., - -^75^5], 



(0f)\- = 5(A:)^'=(A^),, 



% = 4IV (At A,) - 

Gh = -Tr (At A„) + Tr (A^ A.) - -gl 

©E = -10^5^ 



The RGEs for the SUSY breaking parameters in Eq. (|T8D are written as follows 
(47r)^M^(A.)., = {Xu)kAQT)\ + {XuMeT)", + {XuhOH 



+2 



+2 {(A,).,(%)^ + i\,U^T)', + (A.).,%} , 
(4vr)^M^(A.),, = (A.).,(e^)'=, + (A.).,(e^)^,. + (A.),,e^ 

+2 {(A,),,(0^)'=, + iXM^F)'j + {X.^Qh} , 



dM 
d 



+2{{^i),,{eT)'. + i^uerf, + {^U^j: + ©//)} 

+2 {(«:,),,(%)^ + (/€rf).fc(eT)'^ + («:d).,(eE + %)} , 

+2 {(7s:d)fcj(e;p)'=i + (Krf),fc(eT)^j- + {Kd)ij{&j: + e-H 
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d 



(47r)2M 



dM 



2_\i 
N> j 



dM 



+2{\\r[{my)^ + {m^)5i}{\,),, 



where O's are given by 



+2 [4(A:r{K^),' + (m|)4}(Aj);, 
HXir{{ml^)^ + {ml)5i}{X.)i, 

+mr\xi),, + {xir\K),,] 
+f^i{Kr,-2|M5r5i}, 

+2 [5(A:)^'= {(m|^),' + {ml)5l} (A^),, + 5(A:)^'=(A^),, 
26^^771^ + 2 {6Tr (A^m^^A^) + Tr (A|:m|^^A,) 
+IV (A>|^A^) + 3IV (At A„) + IV (At } 

2QTfn]j + 2 {4IV (A^I'^AJ) + 4IV (Atm|^A,) 
+4Tr(A^Aj)} + f,,^(m^-2|M5r), 
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(6^)^. = 2(At)^'=(A,),, + 3(At)^'=(A„)., - -glM,5], 

_ _ _ 94 

(6^)^- = 4(A^)'^AJ)., + (At)^^A.)., - -glM,8], 

= 5(A:)''=(A^),,, 
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=ATi(x\X,)--glM,, (A5d) 
Qh = ^Tt (At A„) + Tr (a^ A,) - -g^M,, (A5e) 
= -lOglM,. (A5f) 



2. Matching conditions at the GUT scale 

In this subsection we show the matching conditions between the SU(5)RN SUSY GUT 
and the MSSMRN at the GUT scale taking account of the dimension five terms in the Kahler 
potential. Although we include only the renormalizable terms in the Kahler potential at the 
Planck scale, higher dimensional terms are induced by the renormalization effects between 
the Planck scale and the GUT scale, because we introduce the higher dimensional operators 
in the superpotential. In order to simplify the treatment, we use logarithmic approximation 
for induced terms in the Kahler potential. For other coupling constants, we explicitly solve 
the RGEs in the previous subsection. 

Up to dimension five terms, the corrections for the Kahler potential in the present model 
are parameterized as follows: 

+kH{H%{m{Hf + kjj{H^ni:)\{H\] + H.C., (A6) 

where we include SUSY breaking parts in the coupling constants using the spurion method 
as follows: 

kx = kx + kxO^ + kxf + kxO^f (X = T, F, H, H) . ( A7) 

At the Planck scale, we assume all the components of Eq. ( [A7| ) are zero. These coupling 
constants are induced from the renormalization between the Planck scale and the GUT scale. 
In the logarithmic approximation, /c's at the GUT scale are given by. 
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^ -2{(A^)M,fc-4(A^),,(«:,),fc- (At)"=(«:,),,.}tG, 
kn ~ -2 {3Tr (aI/s:„) - Tr (At/«.) } to, 
% ^ -2 {Tr (xIko) - 4Tr (A^Ts:,) } tc- 



'G, 



(A8a) 
(A8b) 
(A8c) 
(ASd) 



In the same approximation, kx, kx and kx at the GUT scale are proportional to kx as 
follows: 



kx ~ mo{Ao + AAo)kx, 
kx ~ moAokx, 

kx ~ ml {AliAo + AAo) +2}kx, {X = T, F, H, H) . 



(A9a) 
(A9b) 
(A9c) 

The dimension five terms in Eq. ( |A6| ) modify the normalization of the Kahler potential after 
the SU(5) symmetry breaking. In order to obtain the correct normalization up to O(^), we 
introduce the following new chiral superfields: 



^J.^Y ^ ^j.ryt ^ _1_ ^ eXkrY,} {{^TciTr' - {^TciTT} , (AlOa) 

= {T)a + ^ {(M^- + nhYi] i^TaiFX (AlOb) 

{H'r = {HY + ^{kH + e^kn) i^)\H' (AlOc) 

(H')a = (H)a + l^{kH + 0%) i^)\(Hy. (AlOd) 



Substituting these superfields for Eqs. (|I|) and (^), we can define the coupling constants in 
terms of the new chiral superfields as follows: 



)ij - i'^u)ij - 2{^u)ik{kT)''j - -{K)kj{kTfi + iK)ijkH, 
)ij = + -j{K)ik{kTTj - -{K)kj{kTYi, 



{f^d)ij — {f^d)ij — {^d)ik{kT)^j — {\d)kjik-pY 



il^d)ij — il^d)ij — i^d)ikikT)^j — (XdjijkjY, 



)ij — {^u)ikikp) • {Xu)ijkH- 



(Alia) 
(Allb) 
(Allc) 
(Alld) 
(Alle) 
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The matching conditions for the Yukawa couphng matrices for quarks and leptons are written 
with the above new couphng constants as foUows Q: 

{yu)ij = + ^\^^ii^u')ij + Q{i^u')ijJ , (Al2a) 

{yd)ij = {\d)ij + ^^^Mij - ^(j^d')zj^ , (Al2b) 

iVeh = (Al)i, + ^ [-l{>^d%j - lij^d^,] , (A12c) 

(2/.).. = (A.).,-|m.,. (Al2d) 



Using Eqs. (|^)- (|A12 ) we can relate y's and A's and k's at the GUT scale. The SUSY 



breaking parameters are also defined in terms of the new superfields as follows: 

- ^{\u)^kCkT)''j - ^{^u)kj{kT)\ + {Xu)ijkH, (Al3a) 

iK'k = - \CXu)ik{kT)', + l{\u)kj{kT)^ 

- liXnUkr)', + l{Xu)kAkT)^, (Al3b) 
i^'dh = i^^d)^j - {\d)ik{,kTf^ - (\d)kj{,kp)\ - {Xd)^k{kTfj - {,Xd)kj(kp)\, (A13c) 

^d)ij — 0^d)ij — (^d)ik{kT)'^j — (^d)ijkjY 

- {Xd)^kCkT)'', - {Xd)ijkjr, (A13d) 

{TiDij = - (\y)ik{k-pfj - (\v)ijkH - {\u)ik{k-pfj - {\y)ijkH, (Al3e) 

{k'TY, = CkTY, + {ml)\{kTf,, (Al3f) 

= (M^• + (^F)'fc(M'-, (Al3g) 

k'jj = kH + mj^kH, (Al3h) 

kj{ = k-jj + nijjk-jj. (Al3i) 



^k's in Eqs. ^ of Sec. || should be read as those with prime if we take account of the higher 
dimensional terms radiatively induced in the Kahler potential. 
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The soft SUSY breaking parameters can be expressed using the above new couphng constants 
as follows: 



iyu)^j = + ^ \ ^{K,u')ij + QiK')ij 



+l{XuU4)^ - li)^u)kM)\ + li^uhk}, \ , (A14a) 

2 \^d)ij ~ 2^'^d)i3 



{yd)ij = (Ad)ij + ^ - ^(4) 

-liXdUk^r)', - ^(Ad).,(^^)^ + ^(A,).,%} , (Al4b) 

+ i>^I)kM)\ + ^(Aj)..(fc^)^- + ^(Aj).,%} , (Al4c) 

(y.h = CKh + e {-^iKh + lix.Mk^^)', + liXuWk},^ , (Ai4d) 

= (^tYj + k {{krYj + Ck'^Yj} , (Al4e) 



(m 



2 \ j 

EJi 



(rniY- = [m 



'LJ j 



6 

ml)\-l^{Ck'Ty^ + {ny^}, (Al4f) 

mly^ + ^{CkTy^+{k'^y^}, (Al4g) 

y{CkTy^+i^y^}^ (Ai4h) 



2 _2 



Ir (Al4j) 



ml +ik'H. (Al4k) 



Using Eqs. (|A13| ) and ( |A14|) , we can express the soft SUSY breaking terms of the MSSMRN 



by the input parameters at the Planck scale, Eq. ([T9|) . 



APPENDIX B: APPROXIMATE EXPRESSIONS OF THE PHOTON-PENGUIN 

AMPLITUDES FOR /x ^ 67 PROCESS 

In this appendix we show the explicit forms of the functions which appear in the ap- 
proximated expressions of the photon-penguin amphtudes given in Eq. (P^D- We assume the 
following conditions to derive the expressions: 
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The off-diagonal elements of the slepton mass matrices, {'m'^)^-' and {fn'j^yj are given 
by Eqs. (|25b| ) and (|25d| ), and they are diagonalized with good approximation in the 
basis where and {ycRjij in Eqs. ([2^ ) and ( p3D are diagonal. 



In this basis, the left-right mixing mass of the slepton can be treated as perturbation 
to diagonalize the 6x6 charged slepton mass matrix. 



• The eigenvalues of the slepton mass matrices are almost degenerate and represented 
by rrf'. 

With these conditions, the SUSY contributions to the photon-penguin amplitudes corre- 
spond to Figs. are expressed as Eq. (|53|). In this formula, Oj, and a" are given 
by 



2 

m 



327r 



2 tangly ^(0^)Ai{(07v)vi2 + tan^VK(Oiv)Ai}/2" ( ;;;2- ) \~ 



A=l 



X 



mo(Ae + I AA) + II* tan p ^ 



m 



{'i + \Ao\ytG{tG + tR), 



327r2 cos [3 



j:i0hL)A2{0cR)Air 



2 

m 



m~ 



Xa. 



^)(3+iAon(tG+t 

m 



R) 



(Bla) 
(Bib) 



327r' 



tan^^ ^(0^)^i{(0jv)a2 + tan^^w'(Ojv)Ai}/r 



2 

m 



Xa. 



m J \ m J 5 rriu 



(Blc) 



where 6w is the Weinberg angle. m~o and m~- represent the masses of neutralinos and 

Xa Xa 

charginos, respectively and O^, Ocr and Ocl are unitary matrices which are used to diag- 
onalize the neutralino and chargino mass matrices, M~o and M~_ as follows: 



OnM~oO^ = diag(m~o, m~o, m~o, m~o). 



Ml —mzSw cos P m^s^/sin/? 

M2 mzCvi/cos/3 —mzcw sin P 

-mzSw cos /5 m^cpi/ cos P — /i 

m^Svi/ sin P —mzCw sin /? — 



(B2a) 



(B2b) 
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OcrM-OIji^ = diag(m--,m 



Xl ' X2 



Mz 



M2 \/2miv cos P 
A/2mvi/ sin /? /i 



(B2c) 

(B2d) 
(B2e) 



where rriz and mv^^ are the Z boson mass and the W boson mass, respectively and sw = 
sm9w and cw = cos6'vy. The mass functions /I', and /" in the above formulas are given 
by 



1 

X2 



(1 



{x^ + 4x - 5 + 2(2a; + 1) ln(a;)}, 



{14x3 - 25x2 + I4x - 3 - 2x2 (4x - 1) ln(x)}, 
2x'^(l — x)'* 

^ {5x^ - 37x3 ^ 27x2 + 13x - 



x'^(l — x)' 



6x(7x — 3) ln(x)}. 



(B3a) 
(B3b) 
(B3c) 



APPENDIX C: FCNC EFFECTIVE COUPLINGS IN MSSM 

In this appendix we present the explicit forms of FCNC effective coupling constants for 
j{0 _ ]^o^ j^^ _ j^^ g^Yid Bg — Bg mixings. For — mixing these coupling constants are 
defined in Eq. (|37| ) of Sec. |111 B| and for Bd — Bd and Bs — Bg mixings the coupling constants 
are given by substitution of flavor indices. In the MSSM box diagrams exchanging charged 
Higgs, neutralino, chargino and gluino can contribute to these coupling constants. We first 
define the following neutralino, chargino and gluino vertices for quarks and squarks, 

3 4 6 



^ = E E E {dl {Nf^xPL + Nf^^Pn) x\dx 

+ E E E (cfkxPL + Ct^xPn) XaUx 
■ =1 

+ul {C^^^Pl + C^^xPr) XAdx} 

+ E E E {^fxPL + rf^PR) G^T'^dx 



i=l A=l X=l 



3 2 6 



1=1 A=l X=l 



3 6 8 



1=1 X=l a=l 
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+ui {t^^Pl + rrlPij) c^nx] + H.c, 



(Cl) 



where Pl and Pr are projection operators defined by = (1 — 75)/2 and Pr = (1 + 75)/2 
and is the generator of SU(3) gauge group. The neutrahno-squark couphng constants 
appear in the above formula are given by 



dL 
lAX 



dR 

AX 



jvruL 



uR 
iAX 



-72(72 U tangly (Oat) Ai(f/rf* 



+ 



rridi 



2mw cos (3 



{On)a,{u:\ 



Xi 



-V2g2 



-^(0^)a2 + ^ tan^H/(0^)Ai \ {UDxi 



+ 



rridi 



2m\Y cos P 



(C'iv)^3(?7^)xj+3 



-V2g2 I --tan9wiON)AiiU*)xi+3 



^ muiiVcKu) 



-iO^)A3{U:)xj 



-V2g2 



~{ 2mw sin /5 

lm)A2 + ltanewiO*x)Al]iU:)x^ 



, (^CKm) ^ 

2m^ysin/y 



Xj+3 



(C2a) 



(C2b) 



(C2c) 



(C2d) 



where g2 is the SU(2) gauge couphng constant. Oat is the diagonahzation matrix for neu- 
trahno mass matrix defined in Appendix 0. Ud and Uu are unitary matrices which appear 
in diagonalization of the 6x6 squark mass matrices, m~ and m~ as foUows: 



Udmpi = diag(m| , mj^ , , mj^ , , ) 



mg + mlrrid 

+m| cos2/3(-i + isin^ 61^^)1 



^ cos P{yd + VdfJ'* tanpy 



^ cos P{yd + VdfJ'* tan /3) 



V 



9 t 

mj) + mdm\j^ 



\m\ cos 2(3 siv? 



UumlUl = diag(m?^,m?^,m?3,m2^,m?g,m?J 



(C3a) 
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rriQ + m[m„ 



-^sinp (yu + VufJ'* cot /?) 



V 



sin + cot f3y 
+ |m^ cos 2/5 sin^ 6'^/! 



(C3b) 



where generation indices are suppressed and the mass matrices for down-type and up-type 
quarks are given by {md)ij = mdiS^j and {mu)ij = T^uiiVcKmY j- The chargino-squark cou- 
phng constants in Eq. (|C1|) are given by 



C^AX - 92- 



rridi 



y/2mw cos (3 



j^i v^mi4/sin/5 



iAX 



Ct^X = -92{{OcL)AimX: 
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{OcR)A2{Ud)xj, 



mdi 



:{OcL)A2iU^)xi+3 



(C4a) 
(C4b) 
(C4c) 
(C4d) 



\^mw cos P 

where Ocl and Ocr are the diagonahzation matrices for chargino mass matrix defined in 
Appendix The gluino-squark couphng constants in Eq. ( |C1|) are given by 



rf^ = -V2g,{U:)x^+3, 



iX 



= -V2g,{U:] 



= V2g,{u:: 



Xi+3) 



Xi- 



(C5a) 
(C5b) 
(C5c) 
(C5d) 



where g^ is the SU(3) gauge couphng constant. 

The SUSY contribution to the effective FCNC couphng constants is divided into six parts 
as foUows: 



g = g{H~) + g{H-W) + g{x') + 9{r) + 9{G) + g{Gx') 



(C6) 



where g represents the effective couphng constants in Eq. (^) of Sec. [Iff ij| and its gener- 
alization for — mixing. In the fohowing, coupling constants g are associated with 
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indices — mixing corresponds to i — 1, j — 2 and — {Bg — Bg) mixing 

corresponds to i = 1, j = 3 (i = 2, j = 3). The contribution from box diagrams including 
charged Higgs and up-type quark is given by 

= - W ^ iyCKu)\{VcKM)\{V^KUyk{VcKM)\ 
-•-"^ k,l=l 

X mdjuidi tan^ p d2{m]j- , m]j- , niul, niuf), (C7a) 
X mufmul cof^ P d2{m]j-,m]j-,mul,rnui), (C7b) 

= - W ^ (y^KM)\iyCKM)\{V^KM)\{VcKM)\ 
k,l=l 

X md]muimul dQ{m]j-,m]j-,rnul, m„f), (C7c) 

-■-"^ it,/=l 

X mufmulmdi (io(m|^-,m|-,m4,"^«f), (C7d) 
(i/-) = gL\H-) = 0, (C7e) 



9rl{H ) - -^ 2^ X! (^CKM)M^CKM)'i(VcKM)^fc(^CKM)''i 
k,l=l 



X mdjmufmulmdi do{mH-,m%-,mul, niui), (C7f) 

fc,i=l 

X mdjUiulmdi d2{mH-,m%-,mul, mjf). (C7g) 

where mij- is the charged Higgs mass. The contribution from box diagrams including 
charged Higgs, W boson and up-type quark is given by 

gl{H-W) = 0, (C8a) 



44 



gl{H-W) = iy'cKMyi{VcKM)\{V'cKM)\{VcKM)\ 



V2G_ ^ 

k,l=l 

X rriufmul ' 



cot^ p {d2{mH- , m^, m„^, TUuf) 

- m^(io(m|-,m^,m4,m«D} ' (C8b) 
giniH-W) = = ^lii'(^-W^) = ^fL'(^^-W^) = 0, (C8c) 

gLiH-W) = ^ f (^cW)VVckm)',(K1km)\(Vckm)^ 
xmdjmdirri^ tan^ /3 |(i2(m^- , m^, m„^, m„^) 

2 2 'I 
l^uk'^ul 7 / 2 2 2 2\ I /'rioJ\ 

4m^ — ^o{mH-:rn^:rnuk^mui)> , (C8d) 

- 0. (C8e) 

The contribution from box diagrams including neutralino and down-type squark is given by 

= -198^ ^ ^ ^^^y^^^y 

+NfExN^i^ ^^^^^ do(m^o^ , mjo^ , , ) } , (C9a) 

iZm Up A,B=1 X,Y=l 

X m-^m-o^ do(m|o^,m|o^,m|^,m|^), (C9b) 

/o 4 6 

UStt'Gf ^JBY^^AY [^,BX^^AX ' ^.AX^rBX, 

X m~om~o^ ^o(m|^,m|o^,m|^,m|^), (C9c) 

a/9 4 6 

„S _ ArrfR ATdL*(nTdL prdR* , prdL ATdR*\ 

gRL\X ) — 2/^ 2^ ^^jfiy^^iAy I ^^jBX^^iAX + ^^jAX^^iBX i 

X (i2 (m~o , m~o , m~ ,m~), (C9d) 



giRix') - 
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A,B=1 X,Y=1 



{2iV-,iV-l d,(m|o,m|o 

+Nf^^N^i*^m~o^m~o^ do(m|^, m|o^, m|^)} , (C9e) 

The contribution from box diagrams including chargino and up-type squark is given by 

9lir) = -T:j^j: E C'jl^C^,'5[;Cj^^C;^^irf2K-,m^ m|^,m|J, (ClOa) 
izm Up ji B=i x,Y=i ^ " 

gLix-) = 0, (Clob) 

2 6 

9rR \X ) — ~TOQ^2n 2^ 2_/ '-'JBY'-'iAY'-'jAX'-'iBX 
iZm Up A B=1 X,Y=1 

X mz-mz- dn(m~- , m-- , m- ,m~ ), (ClOc) 
d'RLiD-^^j: E Cf§yCti*Cf^^Cti;, d,{ml m^^^^^ (ClOd) 

D^TT i^F A,B=1 X,Y=1 
/9 2 6 

_s V-c;-^ — V ^ \ " \ ' ^dL r<dL*r<dR r"iR* 

yRL \X ) — ■xnQ^m 2^ 2^ ^JBY^iAY^jAX^iBX 
LZan LfF A,B=1X,Y=1 

X mz-mz- doim-- , m~- , ), (ClOe) 

Xa Xb Xa Xb "X '^Y^ ^ ' 

The contribution from box diagrams including gluino and down-type squark is given by 

^ — 1 

+1m| do(M|,M|,m|^,m|^)}, (Clla) 

^^^(^^ ^ "128^ ]^ , ^rJprf^*rJ^rf|*M| d,{Ml Ml m|^, m|^), (Cllb) 
^^^'^^^ ^ 128^ ^ , ^r^'?rf^*rjf rf|*M| d,(M|, M|, m|^, (Cllc) 



^ T^dRT^dR*T^dL T^dL* T^dLT^dR*T^dRT^dL* 

9rl{'^) - i287r2G^ J^^^ \ \ 3 18 



X d2{MlMlm'~^,m'~;, 
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+^rj^rf^*rjf r^|*M| do{Ml m|, m|^)} , (ciid) 



E 

X,Y=i 



-pdL-pdR*-pdR-pdL* • -pdR-pdR*-pdL -pdL* 

Q^jY^iY ^jx^ix -TgJ-jyJ-iy ^jx^ix 



+^rJ^rf^*rJ«rf|*M| do(M|, M|, m|^, m|^)} , (Clle) 

The contribution from box diagrams including neutralino, gluino and down- type squark is 
given by 

4 6 .9 



1287r2G'i. 



A=l X,Y=i 



I / AjdL -pdL* AidL -pdL* 1 -pdL ]\TdL*T^dL ATdL*\ 



9'rr{Gx' 



9'rr{Gx' 



dx dy' 



(C12a) 



V2 



..on A=l XY=1 ^ ' 



^ ( AT^R T^dL* AjdR T^dL* , j^dR i\TdL*T^dR ATdL*\\ 
"3 y^jAY^ lY ^^jAX^ iX + ^ jY-^^iAY^ jX^^iAX ) J 



dx ' dy ^ 



(C12b) 



V ^ W /\r<iii r<^-f'*r'i^ /\r' 



dL* 
iAX 



I /vrd-R pdL* ^rdi? pdL* 1 pdR AjdL* -pdR AjdL* \ 
^^^jAY^ iY ^^jAX^ iX + ^ jY^^iAY^ jX^^iAX J 

X m-^Ma cio(m|o , M|, m|^), (C12c) 



V I (AjdL T^dR* , T^dL ATdR*\/T^dR ATdL* , Ardi? 

"l287r^Gi7 gl-'^jAy^ iy + ^ jy-'^iAyA^ jx^^iAx + -'^jAx^ ix J 

I / AidR T^dR* AfdL -pdL* i -pdi? Ar(ii2*"pd-£/ -xjdL-^ \ I 
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i^f AjdL -pdL*-pdR ]\TdR* i -potL AjdL* AidR -pdR*\ 



dx' dy' 



(C12d) 



E\ ^ J / /vrdL T^dR* i i-idL a rd-R* \ / /vrdiJ i-idL* i i-idiJ AjdL* \ 
±^on -<^F ]K^^jAY'-iY + jY iAY ) j AX iX 'T^jX^^iAx) 



_L /' AT'^R T^dR* j\jdL -pdL* , -r^dR ]\TdR*T^dL nrdL* \ \ 
+ 2 \^^jAY^ lY ^^jAX^ iX + ^ jY^^iAY^ jX^^iAX ) J 



-L.^ ( Afd'^ T^dL*T^dR AfdR* _i_ pdL nrdL* AidR T^dR*\ 



(C12e) 



The neutralino, chargino, gluino and neutralino-gluino contributions to g^, g'f^ and g^f^' are 
obtained by replacing the sufSx R with L and L with R in the corresponding formulas for 
9rj 9rr 9rr J respectively. The mass functions which appear in the above formulas are 
defined as follows: 



do{x,y,z,w) 



X Iyl{x) 



+ 



y\n{y) 



{y - x){z - x){w - x) {x -y){z -y){w -y) 
^ln(^) w\n{w) 



d2{x,y,z,w) = - 



{x — z){y — z){w — z) {x — w){y — w){z — w) ' 
1 f x^ln{x) ^ y'^ln{y) 



4 [{y — x){z — x){w — x) 
z^ ln(2;) 



+ 



\x - y){z - y){w - y) 
up' \a{yj) I 
{x — z^{y — z){w — z) ' {x — w){y — w){z — w) ] 



+ 



(C13a) 



(C13b) 
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TABLES 
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TABLE L 
Fig. 14. 


Input parameters 


used in the numerical calculation. 


Kib/^cfe and (5i3 are varied in 






LMA 


SMA 



sin^ 26'sun 
sin^ 26latm 



13 

2 f^\r2\ 



Ami2 (eV 
rriu^ (eV) 



1 
1 



1.8 X 10-^ 
3.5 X 10"3 
4.0 X 10^3 



5.5 X 10-3 
1 


5.0 X 10-^ 
3.5 X 10-^ 
2.2 X 10-3 



TABLE IL Parameters for the neutrino sector. LMA (SMA) corresponds to the large (small) 
mixing angle MSW solution for the solar neutrino anomaly. 
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TABLE III. Parameters for Fig. 11 



54 





Bb 


[Bb) 
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S 1 
RL 


Bk {Bk)rl 


{BKfnL 
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TABLE IV. Decay constants and bag 


parameters for — 


and 


— mixing matrix 


elements used in the numerical calculation 
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£k/{£k)sm - 1 


<0.1 


<0.5 






<1 


< 0.05 
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< 0.05 




< 1 


TABLE V. Summary of the SUSY contributions to the observables in 


Fig. 0. 


V" shows 



that some parameter region is excluded by the ^ constraint and hence the branching ratio 
can be just below the present upper bound. means that the SUSY contribution is negligible. 



55 



FIGURES 
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< i i < < i i < 

FIG. 1. Possible large contributions to — > 67 amplitudes, Aj^ and in the present model. 
They are enhanced with a factor mr/mn compared to the other contributions. 
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FIG. 2. A possible large contribution to — >^ 67 amplitude in the present model. 
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3C° f 
FIG. 3. Possible large contributions to — >^ e 7 amplitudes and in the present model. 
The flavor mixing in the left-right mixing is induced by the gauge interaction between the Planck 
scale and the GUT scale through the wavefunction renormalization of 24 Higgs supermultiplet. 

56 



w 



w 




u, c, t 



Sl Sr bR bL 



-im. 





u, c, t 




FIG. 4. One-loop diagrams which contribute to 6 — > s 7 in the SM. The contribution to Cy 
(left) is suppressed by ms/mt, compared to the contribution to C7 (right). 
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FIG. 5. A possible large contribution to 6 ^ 57 amplitude Cy in the SU(5)RN SUSY GUT 
which is not suppressed by m^/mb. 
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FIG. 6. A one-loop diagram which contributes to ek in the SM. A similar diagram which 
contribute to — {Bg — Bg) mixing in the SM can be obtained by replacing quarks in the 
external lines so that s — 6 (s — 6 and d ^ s). 
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FIG. 7. A possible large contribution to er in the SU(5)RN SUSY GUT. There is also a 
crossed diagram because of Majorana nature of gluino. 
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FIG. 8. Possible large contributions to — mixing in the present model. Similar diagrams 
which contribute to Bg — Bg mixing can be obtained by replacing the down quark/squark in the 
diagrams to the strange quark/squark. There are also crossed diagrams. 
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FIG. 9. Branching ratio and P-odd asymmetry of // — > 67 as functions of Be and 9d- 
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FIG. 11. Contour plots of a^^bv^ b(;U ^ e 7), B(t ^ ^ 7), Si^/e^^ - 1, AmB^/Amlf - 1 and 



Am^^/Am^^ — 1 on the ttiq-Mq plane for various choices of the parameters given in Table ID. 
CKM parameters are fixed as 613 = 60° and iFub/Vcbl = 0.08 and fi is taken as positive. Shaded 
regions are excluded experimentally (see text). In (b) and (g), the excluded regions correspond to 
negative stop mass squared are also shown. 
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